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Abstract 

We address known discrepancies between the heavy flavor properties of jets produced at the 
Tevatron collider and the prediction of conventional-QCD simulations. In this study, we entertain 
the possibility that these effects are real and due to new physics. We show that all anomalies can 
be simultaneously fitted by postulating the additional pair production of light bottom squarks 
with a 100% semileptonic branching fraction. 

PACS numbers: 13.85.Qk, 13.20.He, 12.60.Jv 
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I. INTRODUCTION 



In 2001 the CDF experiment reported [1] the observation of an excess of events - con- 
taining a lepton with large transverse momentum, large transverse missing energy, and two 
or three jets one of which (called superjet) contains a displaced secondary vertex and a soft 
lepton due to the decay of presumed b hadrons - over the standard model (SM) expectation 
(17 such events are observed and 4.9 ± 0.6 are expected) 1 . 

The apparent excess of superjets observed by CDF has been modeled in Ref. [3] by postu- 
lating the unconventional production of a low-mass, strongly interacting object that decays 
semileptonically with a branching fraction of the order of 1 and a lifetime of the order of the 
picosecond. Since at the time there were no limits to the existence of a charge —1/3 scalar 
quark lighter than 10 GeV/c 2 [4-6], the supersymmetric partner of the bottom quark was 
a potential candidate. The hypothesis of a very light scalar quark has been investigated in 
many theoretical papers and has prompted careful re-examinations of precision electroweak 
data and e + e~ — > hadrons data [7-21]. It is a fair summary of all this work that no data 
analysis shows convincing evidence for a very light scalar quark. The conjecture itself is not 
favored, but also not completely excluded, in the framework of the minimal supersymmetric 
standard model (MSSM). The CLEO collaboration [22] has searched for pair production of 
light squarks that decay to a lepton, a charmed quark and a scalar neutrino with a 100% 
branching fraction by selecting events that contain two oppositely charged leptons plus a D 
meson below the BB threshold. Reference [22] reports no evidence for the production of 
such squarks with mass between 3.5 and 4.5 GeV/c 2 . However, an exclusion limit depends 
on the modeling of the scalar quark hadronization and decay matrix element, and squarks 
lighter than 3.8 GeV/c 2 could have escaped detection [22]. 

Because of the unresolved experimental situation, the CDF collaboration has performed 
a detailed comparison [23] between measured and predicted heavy-flavor properties of jets 
produced at the Tevatron, intended to search for evidence either supporting or disfavoring 
this conjecture. Reference [23] studies rates of observed and predicted leptons arising from 

1 The DJZ) collaboration has also searched [2] for the presence of such an anomaly, and reports a deficit 
with respect to the SM prediction (2 events observed and 3 ± 0.5 predicted). Unfortunately, because of 
stricter selection requirements and a soft lepton tagging efficiency significantly smaller than that of CDF, 
the resulting 95% C.L. limit to the cross section for producing this type of events is at least a factor of 
two larger than the cross section corresponding to the excess observed by CDF [3]. 
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6-quark decays in jets that recoil against a generic jet or recoil against a jet that also 
contains a lepton. That study finds that in the second case the rate of jets containing a 
lepton from presumed b decays is 50% higher than in the first case. The magnitude of 
the effect is consistent with what is expected in the presence of pair production of light 
squarks with a 100% semileptonic branching ratio. The magnitude of the discrepancy is 
consistent with the result of other measurements performed at the Tevatron. The ratio 
of the \x + fe-jet cross section to the next-to-leading-order (NLO) prediction 2 is measured 
to be 1.5±10% [24], the ratio of the cross section for producing dimuons from presumed 
b- and 6-hadron semileptonic decays to the NLO prediction is 2.6±20% [25, 26], whereas 
this ratio is 1.1±20% when identifying b quarks by locating secondary vertices produced by 
their decays [23, 27, 28]. If one assumes that the above mentioned measurements and their 
quoted uncertainties are correct, the disagreement between data and theory appears to be 
a function of the number of semileptonic decays used to identify b quarks 3 . Reference [29] 
has extended the comparison of Ref. [23] to jets containing a lepton pair. In that study, the 
distributions of the dilepton invariant mass and opening angle are found to be strikingly 
different from those predicted by a conventional-QCD simulation, in which most lepton 
pairs arise from sequential semileptonic decays of single b hadrons. In addition, a study 
of the invariant mass distribution of muon pairs collected by the CDF experiment [30] has 
attempted to rule out 1 bound states of light scalar quarks that are not excluded by the 
SPEAR searches for narrow resonances [5, 31]. Reference [30] improves the SPEAR limit 
by almost an order of magnitude, and sets an average upper credible limit of approximately 
8 eV to the leptonic width of new narrow resonances below the T mesons. An exception is 
the mass of 7.2 GeV/c 2 at which an excess of 250 ± 61 events over a smooth background 
is observed. The statistical significance of the signal (3.5 a) is not enough to establish the 
existence of a new particle. However, the size of this signal is consistent with the production 
of a 1 bound state of spin and charge —1/3 quarks. 

All anomalies listed above certainly require further and independent confirmation. How- 



2 As discussed in Ref. [23], the NLO prediction for producing two central bottom quarks above the same 
transverse momentum threshold has an uncertainty no larger than 15%; in contrast, the NLO prediction 
of the single b cross section has a ~ 50% uncertainty. 

3 This effect could also be a fortuitous result produced by the fact that the distribution of the ratios 
of the bb cross section measurements at the Tevatron, with and without leptons, to NLO prediction has 
approximately a 30% RMS deviation, considerably larger than the quoted experimental uncertainties [28]. 
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ever, we entertain the possibility that all these effects are real, and verify if a new physics 
process can simultaneously model the discrepancies between data and conventional-QCD 
predictions reported in Refs. [23, 29]. Since we use the same data sets, simulated samples, 
and analysis tools we briefly recall those studies. Reference [23] compares a QCD simulation 
based upon the herwig [32] and QQ [33] Monte Carlo generator programs to a data sample, 
referred to as inclusive lepton sample, that consists of events with two or more jets, one of 
which is required to contain a lepton with pt > 8 GeV/c (trigger lepton); this jet is referred 
to as lepton-jet, whereas jets recoiling against it are called away-jets. That study uses mea- 
sured rates of lepton- and away-jets containing displaced secondary vertices (SECVTX tags) 
or tracks with large impact parameter [34] (JPB tags) in order to determine the bottom and 
charmed content of the data. In order to remove the uncertainty of the predicted heavy 
flavor production cross sections, the simulation is tuned to match the heavy-flavor content 
of the data. In Ref. [23] rates of SECVTX and JPB tags, as well as the event kinematics, are 
well modeled after tuning the parton- level cross sections predicted by the herwig generator 
within the experimental and theoretical uncertainties. However, the data have a 50% excess 
of events in which one additional jet (away-jet) contains a soft lepton with pt > 2 GeV/c 
(SLT tag). Reference [29] has extended the comparison between data and tuned simulation 
to lepton-jets that, in addition to the trigger lepton, contain an additional soft lepton with 
opposite sign charge (OS). After removing the background arising from fake leptons by sub- 
tracting rates of lepton pairs with same sign charge (SS), the yields of OS-SS lepton pairs 
are compared to the tuned simulation. For invariant masses smaller than 2 GeV/c 2 , the 
data are found to be largely underestimated by the simulation. 

In the present study we try to model the same data used in Refs. [23, 29] by implementing 
the conventiona-QCD simulation, referred to as SM simulation in the following, with the hy- 
pothetical pair production of sbottom quarks ( b s ). Section II describes the implementation 
of this new process into the SM simulation. In Sec. Ill we tune the SM+6 S simulation to also 
reproduce measured yields of jets with SLT and dilepton tags. Section IV investigates the 
sensitivity of the data to other models. Section V presents a crude estimate of the lifetime 
of the hypothetical object causing the excess of SLT tags. Our conclusions are presented in 
Sec. VI. 
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II. PRODUCTION AND DECAY OF SBOTTOM QUARKS 



The production of b s quarks that decay to I X with a 100% branching ratio would be 
sufficient to model the excess of away-jets with SLT tags reported in Ref. [23]. In order to 
model the excess of jets containing a lepton pair reported in Ref. [29], we need an extra 
source of leptons. This could be provided by the sequential semileptonic decay of c quarks, 
and we choose the decay b s — > c I u s , where v s is a massless non-interacting scalar. 

We model the b s decay with the matrix element discussed in more detail in Appendix A. 
In this model, the decay is mediated by the higgsino coupling to the right-handed matter. 
This matrix element produces lepton pairs from sequential decay that have an invariant 
mass distribution different from that produced by the V — A matrix element of conventional 
heavy-flavor decays (see Fig. 1). We implement the production of pairs of scalar quarks in 
pp collisions into the herwig generator with the method described in Appendix B. The 
method starts with the generation of bb pairs with a quark mass of 3.65 GeV/c 2 and a 
lifetime of 1 ps. At the end of the 6-hadronization process, the formed B hadrons are turned 
into fictitious B hadrons. We have generated samples of simulated events with the same 
luminosity of the conventional-QCD samples used in Refs. [23, 29]. The NLO calculation of 
the process pp — > b s b s , implemented into the prospino Monte Carlo generator [35], predicts 
a cross section that is approximately 15% of the NLO cross section for producing pairs of 
quarks with the same mass [36]. Therefore the parton- level b s cross section in the HERWIG 
simulated samples is overestimated by a factor of seven. 
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FIG. 1: Distributions of the invariant mass M and the opening angle 9 between leptons arising 
from the sequential decays b — > I c with c — > Z s f and 6 S — > Z c i/ s , with c — > Z s za We set 
m c = 1.5 GeV/c 2 . Leptons are not simulated through the detector. 
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III. FIT OF THE DATA WITH THE SM SIMULATION IMPLEMENTED WITH 
PAIR PRODUCTION OF SCALAR QUARKS 

In Ref. [23], after tuning the SM simulation to fit the observed rates of lepton- and 
away-jets with SECVTX and JPB tags due to heavy flavor, we have observed a discrepancy 
between measured and predicted rates of away-jets with soft lepton tags. In this study, 
we use a simulation that includes pair production of b s quarks with 100% semileptonic 
branching ratio; we tune this production by also fitting rates of away-jets with SLT tags. 
Rates of tagged jets in the data are taken from Tables II and VII of Ref. [23]. We also 
fit rates of (OS-SS) lepton pairs (Dil) contained in the same jet, taken from Table III of 
Ref. [29]. We use the same notations and fitting procedure described in Sec. VIII of Ref. [23]. 

In the fit we tune the parton-level cross section of the different heavy-flavor production 
mechanisms [direct production (LO), flavor excitation, and gluon splitting] to reproduce 
the rates of tags observed in the data. Following the procedure of Ref. [23], we use two 
fit parameters to account for the uncertainty of the luminosity and selection cuts of the 
electron and muon data; these factors include the uncertainty of the b direct production. 
Five fit parameters are used to determine the relative weight of the remaining parton-level 
cross sections for producing conventional heavy flavors (b and c) with respect to the b direct 
production. These five parameters are constrained to the default values as described in 
Sec. VIII of Ref. [23]. Direct production and flavor excitation of b s quarks are weighted 
with two additional fit parameters. As explained in Appendix B, the value of the flavor 
excitation cross section returned by the fit also includes gluon splitting contributions. The 
efficiency for finding SECVTX and JPB tags in the simulation have been tuned using a 
subset of the inclusive lepton sample [37]. Since we are now entertaining the hypothesis that 
this data sample might be contaminated by unconventional production, in the fit we also 
weight the b- and c-quark tagging efficiencies with free parameters (scale factors) to allow 
for the possibility that they might be wrong. We assume that the simulated efficiency for 
tagging b s hadrons is correct 4 . 

The result returned by the best fit with the SM+6 S simulation is listed in Table I. This 

4 We have arbitrarily fixed the b s lifetime; the introduction of a scale factor for the simulated b s tagging 
efficiency is equivalent to introducing an uncertainty on the lifetime. The dependence on the lifetime is 
studied in Sec. V. 
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table should be compared to Table IV of Ref . [23] . From the fit parameters, we derive a cross 
section for producing b s quarks that is (11 ± 2)% of the production cross section of b quarks 
with the same mass, in reasonable agreement with the theoretical prediction mentioned in 
the previous section. Table II compares tagging rates in the data and in the simulation 
tuned according to the best fit. This table should be compared to Tables VI and IX of 
Ref. [23] and Table III of Ref. [29]. 

Because of the approximations in the modeling of the b s production and decay listed in 
Appendix B, a determination of the systematic errors of the b s production cross section is 
beyond the scope and reach of this study. As an example, the herwig spectator model 
predicts a reasonable, but also quite arbitrary, B mass of approximately 4 GeV/c 2 . A 
±0.3 GeV/c 2 change of the B mass results in a ±40% change of the b s production cross 
section returned by the best fit to the data. 

Figures 2 to 4 compare observed distributions of the invariant mass and the opening angle 
of dileptons in the same jet to the SM±6 S simulation. The simulation is normalized with the 
fit parameters listed in Table I. These data appeared to be anomalous when compared to the 
SM simulation in Ref. [29]. For completeness, Figure 5 compares the transverse momentum 
distribution of dileptons in the same jet. In conclusion, the hypothesis of a light parton with 
a 100% semileptonic branching ratio and a production cross section typical of scalar quarks 
fits the yet unresolved discrepancies between the heavy flavor production observed at the 
Tevatron and the prediction of conventional-QCD simulations reported in Refs. [23-26, 29]. 



IV. SENSITIVITY OF THE DATA TO OTHER MODELS 

In the previous section we have used a b s mass of 3.65 GeV/c 2 . This value was chosen 
because, in the invariant mass spectrum of muon pairs collected by the CDF experiment, 
we have seen an excess of events at a mass of 7.22 GeV/c 2 consistent with the production of 
a bound state of scalar quarks. Since this excess could be a statistical fluctuation, we have 
also generated several simulations with b s masses ranging from 2.5 to 5.0 GeV/c 2 . We find 
that the data are reasonably well modeled when using any b s mass between 3 and 4 GeV/c 2 . 

As correctly pointed out in Ref. [9], the measured semileptonic branching fraction of b 
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TABLE I: Result of the fit with the SM simulation implemented with the pair production of scalar 
quarks b s . The scalar quark mass is set to 3.65 GeV/c 2 and its lifetime to 1.0 ps. The fit yields a 
X 2 of 10 for 20 degrees of freedom. 



SECVTX scale factor (b quarks) 


0.98 ± 0.03 


SECVTX scale factor (c quarks) 


1.02 ±0.22 


JPB scale factor (6 and c) 


1.00 ±0.02 


e norm. 


0.96 ± 0.06 


/i norm. 


u.yo it u.ud 


c dir. prod. 


1.01 ±0.14 


b flav. exc. 


0.79 ±0.16 


c flav. exc. 


0.87 ±0.26 


g^bb 


1.32 ±0.18 


9 -> cc 


1.36 ±0.34 


b s dir. prod. 


0.07 ±0.03 


b s flav. exc. + GSP 


0.20 ± 0.04 



hadrons would be higher than the standard model expectation if m ba + m Us < m b 5 . We 
have produced simulations of the b s decay with v s masses of the order of 1 GeV/c 2 . The 
predicted distributions for the dilepton kinematics look very similar to the ones shown in 
Figs. 2 to 4; the fit to the observed tagging rates yields a higher b s production cross section 
(17% of that for producing b quarks with the same mass). 

We have chosen a b s semileptonic decay mediated by the right-handed component of the 
Higgsino coupling. We have also explored the use of a decay matrix mediated by the chargino 
left-handed coupling to ordinary matter (see Appendix A) and the hypothesis of a -R-parity 
violating b s — > I c decay. Both models produce leptons with a stiffer energy spectrum than 
the right-handed decay mode, but also provide a fair description of the data when using 
b s masses approximately equal or slightly smaller than 3.2 GeV/c 2 . In these cases, the b s 

5 In this case the decay b — > b s Dv a with b s — > IX is allowed and does not suffer CKM suppression. This decay 
would however produce leptons with momentum significantly softer than those from direct semileptonic 
b decays. 
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TABLE II: Numbers of tags due to heavy flavors in the data and in the SM+6 S simulation. The 
simulation is normalized according to the fit result listed in Table I. The SM (6 and c quarks) 
and b s contributions are also listed separately. The electron (muon) sample consists of events in 
which the trigger lepton is an electron (muon). In the data, a fraction of the trigger leptons is 
faked by hadrons that mimic the lepton signature: the heavy flavor purity determined by the fit is 



(44.5 ± 1.3)% for the electron 


sample and (58.5 


! ± 2.2)% for the 


muon sample. 










Electron sample 






Tag type 


Data 


Simulation 


SM 


6 s -dir 


6 s -(f.exc + GSP) 


HFl-jet 


68544 


30478.0 ± 878.4 


26040.0 ± 1282.2 


575.7 ± 184.9 


3862.4 ± 882.7 




73335 


32559.0 ± 943.7 


27942.2 ± 1354.9 


593.1 ± 190.5 


4023.7 ± 919.5 


HFT, SEC f 

I jet 


10115.3 ± 101.7 


10135.9 ± 144.7 


8641.1 ±328.2 


189.0 ±61.4 


1305.8 ± 302.6 




11165.4 ± 115.8 


11136.1 ± 142.1 


9279.5 ± 379.7 


235.0 ± 75.9 


1621.7 ± 372.1 


HFTSEC t 


3719.8 ± 93.5 


3678.2 ± 102.8 


3384.1 ± 124.4 


173.8 ± 56.5 


120.3 ± 28.8 


HFTJPf et 


4021.3 ±140.7 


3972.5 ± 104.1 


3572.0 ± 130.4 


228.2 ± 73.7 


172.2 ±40.4 


HFDT SEC 


1375.2 ± 37.6 


1373.4 ±55.9 


1279.0 ± 59.5 


53.9 ± 18.0 


40.5 ± 10.7 


HFDT JPB 


1627.8 ±43.7 


1647.2 ± 53.6 


1458.9 ± 57.8 


89.7 ± 29.4 


98.6 ± 23.7 


HFT a-jet 


862.2 ± 114.8 


840.8 ±98.1 


534.7 ± 63.0 


172.0 ±58.0 


134.1 ± 34.3 


fjprpSLTSEC 


322.1 ± 23.3 


309.6 ±26.9 


219.8 ±20.7 


53.1 ± 17.7 


36.6 ±9.6 


jjprpSLTJPB 


349.4 ± 26.3 


353.0 ±29.9 


224.5 ±20.1 


69.2 ±22.9 


59.3 ± 14.8 


Dil 


1111 ± 54.2 


1149.6 ± 130.7 


752.2 ±88.0 


52.2 ± 17.8 


345.2 ±86.9 


Dil SEC 


401.0 ± 25.0 


436.7 ± 35.9 


312.5 ± 29.1 


14.2 ±4.9 


109.9 ±27.1 


Dil JPB 


464.0 ± 26.8 


480.6 ±38.9 


323.1 ± 29.0 
Muon sample 


18.1 ± 6.1 


139.4 ±33.7 


Tag type 


Data 


Simulation 


SM 


fe s -dir 


6 s -(f.exc + GSP) 


HFl_j et 


14966 


o i yu . u zn o ±u. o 


71 fil AOZ 7 

/ 1U1.1 ZL\ *±£0' I 




J.'iU-L.O 011.1 


HF a -j et 


16460 


9557.2 ± 335.5 


7831.5 ±456.7 


246.5 ± 79.0 


1479.1 ± 328.4 


HFT SEC 
I— jet 


3657.3 ± 60.8 


3635.2 ± 81.2 


2962.3 ± 137.5 


94.9 ±30.8 


578.0 ± 130.2 


l— jet 


4068.6 ± 66.2 


4051.1 ± 83.6 


3234.5 ± 156.8 


115.3 ±37.2 


701.3 ± 156.3 


HFT SEC 


943.8 ± 35.2 


957.2 ±38.7 


842.3 ±43.5 


74.5 ± 24.2 


40.4 ±9.7 


HFT .JPB 


1086.8 ±45.0 


1052.0 ±40.0 


899.2 ±45.2 


95.8 ±30.9 


56.9 ± 13.2 


JjrpSEC 


452.6 ± 21.6 


467.4 ± 27.6 


416.6 ±28.5 


30.3 ± 10.1 


20.5 ±5.3 


DT JPB 


546.4 ± 25.1 


557.3 ±25.5 


468.8 ± 26.9 


46.5 ± 15.2 


42.1 ± 10.0 


HFT SLT 

a— jet 


271.9 ± 34.9 


235.3 ±31.3 


127.6 ± 18.2 


67.8 ±22.9 


39.8 ± 10.2 


uprpSLT-SEC 
nr 1 a-jet 


63.2 ± 10.0 


82.3 ± 10.4 


46.8 ± 7.5 


23.8 ± 7.9 


11.7 ± 3.2 


TT iprpSLT JPB 


103.0 ± 12.4 


101.9 ± 11.8 


54.9 ± 7.7 


29.9 ±9.9 


17.0 ±4.3 


Dil 


336.0 ± 28.5 


371.1 ± 45.9 


232.9 ±31.6 


23.1 ± 7.9 


120.8 ±29.9 


Dil SEC 


158.2 ± 15.9 


168.3 ± 17.2 


116.9 ± 15.0 


8.0 ± 2.7 


43.4 ± 10.6 


Dil JPB 


171.5 ± 16.5 


177.7 ± 17.1 


115.5 ± 13.9 


9.5 ±3.2 


52.7 ± 12.6 
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FIG. 2: Observed distributions of the invariant mass, M, and the opening angle, 9, of dileptons 
in the same jet are compared to the simulation without (SM) and with 6 s -quark production. The 
simulation is normalized with the fit parameters listed in Table I. 

production cross section returned by the best fit to the observed tagging rates is ~ 5% of 
the production cross section of quarks with the same mass. We note that a sbottom quark 
with a mass of approximately 3.2 GeV/c 2 is not favored by the study in Ref. [30], but is out 
of the mass range investigated in the CLEO analysis [22]. In the next two subsections, we 
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FIG. 3: Observed distributions of the invariant mass M and the opening angle 8 of lepton pairs 
contained in a jet tagged by the SECVTX algorithm are compared to the simulation. 



present two additional tests of our model. 
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FIG. 4: Observed distributions of the invariant mass M and the opening angle 8 of lepton pairs 
contained in a jet tagged by the JPB algorithm are compared to the simulation. 

A. Study of the l~ D° system 

In Sec. X B of Ref. [23], the b purity of the inclusive lepton sample has been verified by 
counting the number of D° mesons accompanying the trigger leptons. In the simulation, 
most of the l~ D° (and its charge conjugate) pairs arise from decays of single b hadrons. 
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FIG. 5: Transverse momentum distributions of the trigger (lep) and soft (SLT) leptons contained 
in the same jet are compared to the simulation. 



In that study, the ratio of the number of identified l~ D° pairs in the data to that in the 
SM simulation is 0.92 ± 0.13. This agreement has favored the initial assumption that b s 
quarks decay into c quarks because our model replaces a fraction of the bb production with 
b s b s production. We search the SM+6 S simulation for D° — > K~n + candidate around the 
direction of the trigger lepton following the same procedure described in Sec. X B of Ref. [23] . 



14 



In the mass range of 1.82 — 1.92 GeV/c 2 we identify 147 ± 15 D° mesons (after subtracting 
a background of 47 ± 6 events), in agreement with the number 126 ± 15 found in the data 
(after removing a background of 79 ± 6 events). In the SM+6 S simulation 24% of the l~ D° 
pairs arise from b s decays. Since the b s quark is assumed to be lighter than a b quark, l~ D° 
pairs from b s decays should cluster at smaller invariant masses than those from b decays. 
It is therefore of interest to compare the l~ D° invariant mass distribution in the data and 
the simulation (see Fig. 6). In this comparison the shape of the background is evaluated 
by using wrong-sign combinations {l~ K + n~) with the K + n~ invariant mass in the range 
1.82 — 1.92 GeV/c 2 . Figure 7(a) compares l~ D° invariant mass distributions in the data 
and in the simulation after background removal 6 . 

A comparison of the two distributions in the mass range 2.4 — 5.2 GeV/c 2 yields a 
X 2 of 10 for 13 DOF. However, the SM part of the simulation contains 112 ± 14 l~ D° 
candidates, which is also in agreement with the data (126 ± 15 candidates). The invariant 
mass distribution of the l~ D° candidates contributed by the SM part of the simulation, 
shown in Figure 7(b), also agrees with the data (the x 2 of the comparison of the invariant 
mass distributions is 6 for 13 DOF). In the limited statistics of this sample, the l~ D° data 
do not provide additional support to our model. 



6 Because of the pr cuts used to select leptons, pions, and kaons, the pr threshold of the I D° system 
decreases rapidly with increasing I D° invariant mass. Therefore, the l~ D° invariant mass distribution 
is also affected by the pr distribution of the B mesons that in turn is affected by the shortcomings of the 
herwig fragmentation and hadronization model (see Appendix B). 
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FIG. 6: Invariant mass distributions of Z~ D° candidates in (a) the data and (b) the SM+6 S simu- 
lation. The dashed histograms represent the background evaluated using wrong-sign combinations 
(see text). 
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TABLE III: Numbers of J/ip mesons in the data and in the SM+6 S simulation. 



J/1> 



e ' e 



DU J / B 



data 
176.0 ± 14.4 
57.8 ± 8.8 
61.2 ±8.4 



simulation 
131.4 ±18.7 
61.3 ±9.5 
57.5 ±8.3 



J/ip -» fi+fi 
data simulation 
83.0 ±9.4 48.3 ±8.6 

31.9 ± 5.8 25.7 ±5.2 

29.6 ±5.7 28.1 ±5.2 



B. Rates of lepton-jets containing J/ip mesons 

J/ip mesons are produced in the decay of B hadrons but not in the decay of B hadrons 
with a 100% semileptonic branching ratio. The SM±6 S simulation predicts a rate of J/ip 
mesons from B decays that is approximately 20% smaller than that predicted by the SM 
simulation (see Table XIV of Ref. [23]). Rates of J/ip in the data and in the SM+6 S simu- 
lation are listed in Table III. Rates of J/ip mesons contained in lepton-jets with SECVTX 
and JPB tags arise only from B decays, and are well predicted by the SM±5 S simulation. 
Before tagging, (50 ±5)% of the J ftp mesons identified in the data (first row of Table III) 
do not arise from B decays [23, 38]. The SM±5 S prediction (180±20) is still slightly larger 
than the data in which 130 ±11 J/ip mesons originate from B decays. 

V. LIFETIME 

In this section we use several techniques to estimate the lifetime of the object producing 
the excess of soft lepton tags with respect to the SM simulation. We extract the lifetime 
by a data to simulation comparison of the distributions of the following quantities strictly 
related to the lifetime: 

1- r dii = Lxy — , where M and pt are the invariant mass and transverse momentum of 
a lepton pair (trigger lepton plus soft lepton), and L xy is the projection on the pt 
direction of the distance between the primary vertex of the event and the vertex 
determined by the lepton tracks. 

A A ~~* — * 

2. T\ — 7re |^| , where d is the vector corresponding to the distance of closest approach of 
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FIG. 7: The distribution of the I D invariant mass after background removal is compared to the 
prediction of (a) the SM+6 S simulation and (b) the SM contribution only. 
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a lepton track to the primary event vertex, and pr is the transverse momentum of the 
jet containing a trigger (lep) or soft (SLT) lepton. 

We test these methods using lepton pairs consistent with J/ip mesons produced by b decays 
(Dil^ PB in Table III). Figure 8 shows agreement between data and simulation within the 
limited statistics of the J/ip sample. 




FIG. 8: Test of the distributions used to estimate the b s lifetime using J/ip mesons contained in a 
jet tagged by the JPB algorithm (see text). 

Figures 9 and 10 compare lifetime distributions in the data and in the simulation imple- 
mented with the production of a scalar quark with a 1.0 ps lifetime and normalized as in 
Table I. Using the track parameters of the trigger and soft leptons in the same jet we con- 
struct the quantities r^i and n ep . Using soft leptons in away-jets with JPB tags, we construct 
the quantity r SLT . In the data this last distribution has been obtained after subtracting the 
expected amount of fake SLT tags. The t S lt distribution for fake tags is approximated using 
the distribution of all SLT candidate tracks in these jets and is shown in Fig. 11. For each 
distribution we evaluate the quantity 

X Y ED(i) 2 + ES{i) 2 
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FIG. 9: Lifetime distributions derived using lepton pairs in the same jet are compared to the 
prediction of the simulations without (SM) and with b s production. The simulation is normalized 
according to the fit result listed in Table I. A comparison of the measured and predicted T\ cp dis- 
tributions yields a x 2 of 8 for 15 DOF (—1.0 < T\ cp < 2 ps); the comparison of the distributions 
yields a x 2 of 8 for 12 DOF (-0.4 < r di i < 2 ps). 
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FIG. 10: Lifetime distributions derived using soft leptons contained in away-jets with JPB tags. 
The comparison of the tslt distributions in the data and the SM+6 S simulation yields a x 2 °f 15 
for 10 DOF (0 < tslt < 2 ps). 

where D(i) and S(i) are the contents of the i th -bin of the measured and predicted distribu- 
tions, respectively (ED(i) and ES(i) are the corresponding errors). The x 2 is constructed 
avoiding large negative and positive tails which are not well modeled by our detector simu- 
lation [37]. We have produced a set of scalar quark simulations with lifetime values ranging 
from 0.2 to 3.0 ps which we add to the SM simulation. We first fit each of these simulations 
to the data following the same procedure used in Sec. II. Since some tagging rates depend 
upon the lifetime, we use the x 2 value returned by the fit as an additional discriminant. The 
X 2 yields as a function of the b s lifetime are shown in Fig. 12. 

The fit of the measured tagging rate yields a lifetime of 0.54 ± 0.3 ps 

The Ti e p distribution yields a lifetime of 1.0 ± 0.5 ps 

The Tdii distribution yields a lifetime of 0.9 ± 0.6 ps 

The tslt distribution yields a lifetime of 1.55^o;4 ps 



21 




FIG. 11: The lifetime distribution of soft leptons in away-jets tagged by JPB is compared to the 
expected contribution of fake SLT tags. 



The average of the four results is 1.0 ps; the uncertainty defined as the RMS deviation 
of the four measurements is 0.4 ps. 



VI. CONCLUSIONS 



We review a number of known discrepancies [23-26, 29] between the heavy-flavor proper- 
ties of jets produced at the Tevatron and the prediction of conventional-QCD simulations. 
These discrepancies may well be the result of larger than estimated experimental difficulties 
hidden in all bb cross section measurements [28]. However, in this study, we entertain the 
possibility that these effects are real, and investigate the constraints imposed by the data on 
the hypothesis that a fraction of the heavy flavor hadrons produced at the Tevatron might 
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FIG. 12: Yields of x 2 resulting from comparisons of data and simulation as a function of the b s 
lifetime. 



be due to unconventional sources. We postulate the existence of a sbottom quark b s that 
decays to a lepton, a charmed quark, and a neutral scalar non-interacting particle v s with 
a 100% semileptonic branching fraction, and implement its production into a conventional- 
QCD simulation based on the herwig Monte Carlo generator program. We initially use a b s 
mass of 3.6 GeV/c 2 since it is not excluded in the study of Ref. [30], a massless u s , a lifetime 
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of 1 ps, and a decay mediated by the higgsino coupling to the right-handed matter. We 
use this generator to simulate the inclusive lepton sample collected by the CDF experiment 
in the 1992 — 1995 collider run. The data sample consists of events with two or more jets, 
one of which is required to contain a lepton with > 8 GeV/c; this jet is referred to as 
lepton-jet, whereas jets recoiling against it are called away-jets. Following Ref. [23], this 
study uses rates of lepton- and away-jets containing displaced secondary vertices (SECVTX 
tags), tracks with large impact parameters (JPB tags), and soft (pr > 2 GeV/c) leptons 
(SLT tags) in order to determine the bottom, charmed, and b s content of the data. We 
tune the parton-level cross section predicted by the simulation within the experimental and 
theoretical uncertainties. The best fit to the data returns a b s production cross section at 
the Tevatron that agrees with the NLO prediction for pair production of scalar quarks. Ac- 
cording to the fit, the b s contribution to the inclusive lepton sample is approximately 23% 
(35%) of that of the b quark in the electron (muon) sample. This contribution models the 
invariant mass and opening angle distribution of lepton pairs contained in the same jet that 
were found anomalous when compared to a SM simulation in Ref. [29] . We have studied the 
sensitivity of the data to slightly different models and investigated the lifetime of the hypo- 
thetical b s quark. The best fit to the data replaces approximately 20% of the b production 
with b s production. In our hypothetical model, the invariant mass distribution of I D° pairs 
coming from b s decays is expected to peak at smaller values than that of I D° pairs coming 
from b decays. In the limited statistics of the data, we see no clear evidence of this. Since b s 
decays do not produce J/ip mesons, we have verified that rates of J/ip mesons in the data 
are consistent with those predicted by the simulation that includes b s production. 



VII. ACKNOWLEDGMENTS 



We thank the Fermilab staff, the CDF collaboration, and the technical staff of its par- 
ticipating Institutions for their contributions. We thank M. Spira for making available the 
code of the prospino Monte Carlo generator, and M. Mangano for working out many de- 
tails of the matrix element used to decay sbottom quarks. This work was supported by the 
U.S. Department of Energy and the Istituto Nazionale di Fisica Nucleare. 



24 



APPENDIX A: DECAY b s -> \cu s 

The decay of 6 S hadrons is modeled using the spectator model (routine hwdhvy) available 
in the version 5.6 of herwig. In this model, herwig normally weights the phase space of 
the three-body semileptonic decay of a b quark with the V-A matrix element. In the case 
of a scalar quark we have modified the HERWIG code and weight the phase-space with the 
matrix element: 

K[(l - z c )(l - Zl ) - R Us + R c (z c - Zl + R„ s - R c )\ (Al) 



dz c dzi 

where K is a normalization factor, R c = m 2 c /ml s1 and R Us = m% /m% . 
Here z c and zi are defined as 

2Pb s ■ Pc , 2pb s ■ Pi 
z c = ^ — and z\ = 5 — 

m b s m b s 

The phase space limits are: 

2^R C <z c <l + R c -R Us 

1 + R c — R u , — z c 1 + R c — R Us — z c 

< zi< 



l-[z c - yjzl - AR c }/2 l-[z c + ^Jz 2 c - 4R c }/2 

This matrix element for the decay b s — > clv a is derived from the tree level calculation outlined 
in Ref. 



M 



-g 2 v 3 ^\v ci \ 



(Pi - Pc) 2 
j=i 



m b Uj 2 sind 

Uj\cos6 



\/2m\Y cos /3 



m c V* 2 cos9 



m Xi u (Pc)-Pflf(p;) + -= —u(Pc)p / dPrv(pi) 

V ifnw sin/3 



(A2) 



where x/ s are the chargino mass eigenstates, and U and V* are the mixing matrices for 
the right and left-handed charginos, respectively. The first subscript of U, V corresponds 
to mass eigenstates and the second to weak eigenstates (1 for the gaugino and 2 for the 
higgsino). Here tan/3 is the ratio of the vacuum expectation values of the two Higgs fields, 9 
is the mixing angle between left-handed and right-handed scalar quarks, F L is the fraction 
of left-handed component of the scalar neutrino, and |V^| is the CKM matrix element. 
Equation (Al) follows from (A2) when the decay is mediated by the higgsino coupling to 
the right-handed matter. If the decay is mediated by the gaugino coupling to the left-handed 
matter the matrix element is 

K{z c + Zl -l + R Us - R c ). (A3) 



dz c dz\ 
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In the latter case the two fermions in the final state are both left-handed and tend to 
be produced back-to-back since the initial state is spinless. Using equation (Al) the two 
fermions in the final state have opposite handedness and tend to be produced more collinear 
than when using equation (A3). It follows that the matrix element (Al) produces leptons 
with a momentum distribution appreciably softer than that produced by the matrix element 
(A3) or by a phase space decay; this latter mode has been used in Ref. [22] to search evidence 
for such a scalar quark. 

The c quark emerging from the b s decay is recombined with the spectator quark by the 
herwig spectator model. The decay of the resulting c hadron is performed with the QQ 
Monte Carlo program. In the spectator model, excited D-meson states are produced only 
by the hadronic current carried by the virtual W (so called upper vertex). Since we impose 
that the gauge fermion involved in the scalar quark decay has only leptonic decay modes, 
the simulation produces a very simple list of c hadrons with respect to the QQ generator (see 
Fig. 13). The value of the B masses returned by the spectator model, which does not have a 
look-up table similar to that for B hadrons, are no more than an educated guess, probably 
with an uncertainty of a few hundred MeV. The fractions of different c-hadrons produced 
by the spectator model are also quite arbitrary. 
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FIG. 13: Mass spectrum of: b and 6 S hadrons formed by the herwig generator using nib = 
4.75 GeV/c 2 and = 3.65 GeV/c 2 ; c hadrons from 6-hadron decays modeled with the QQ 
generator; and c hadrons from 6 s -hadron decays modeled with the herwig spectator model. 



APPENDIX B: MODELING OF THE b s PRODUCTION 

We implement the production of pairs of scalar quarks in pp collisions using option 1705 
of the version 5.6 of the herwig generator program. This option generates bb events with 
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LO and flavor excitation diagrams, but without gluon splitting diagrams. We use a 6-mass 
value of 3.65 GeV/c 2 and initially set the b s lifetime equal to 1.0 ps. The B hadrons, formed 
at the end of the 6-fragmentation process, are changed into fictitious B hadrons and decayed 
by the herwig spectator model fed with the b s — > I c v s matrix element described in the 
previous Appendix. The decay of the c-hadron formed by herwig using the spectator quark 
is modeled with the QQ generator program [33]. Our modeling of the b s production suffers 
of the following approximations: 

1. The b s b s production cross section has not the same dependence on the velocity (3 of 
the partons in the final state as the bb cross section [40]. However, we study jets 
with transverse energies larger than 15 GeV and this alleviates the problem, since the 
corresponding partons have energies larger than 20 GeV. 

2. The spectrum of gluon emission in the hadronization process of a spin quark differs 
from that of a spin 1/2 quark. However, the average energy loss due to perturbative 
gluon emission off a spin-0 and a spin- 1/2 particle is believed to be small [41]. 

3. Option 1705 does not evaluate the gluon splitting contribution to b s production. How- 
ever, we tune the b s production to reproduce the observed rates of tagged lepton- and 
away-jets. When only looking at yields of tagged lepton- and away-jets, Table X of 
Ref. [23] shows that events arising from flavor excitation and gluon splitting are quite 
similar. Therefore we fit the data with the SM simulation implemented with the scalar 
quark production leaving its direct production and flavor excitation cross sections as 
additional free fit parameters; in first approximation, the flavor excitation cross section 
determined by the fit also accounts for the gluon splitting contribution. 
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